Protostars, Dust Globules, and a Herbig-Haro Object in the LMC 

Superbubble N 51D 



Y.-H. Chu 1 , R. A. Gruendl 1 , C.-H. R. Chen 1 , B. A. Whitney 2 , K. D. Gordon 3 , L. W. 
Looney 1 , G. C. Clayton 4 , J. R. Dickel 1 , B. C. Dunne 1 , S. D. Points 5 , R. C. Smith 5 , and 

R. M. Williams 1 

ABSTRACT 

Using Spitzer Space Telescope and Hubble Space Telescope observations of the 
superbubble N51D, we have identified three young stellar objects (YSOs) in dust 
globules, and made the first detection of a Herbig-Haro object outside the Galaxy. 
The spectral energy distributions of these YSOs suggest young massive stars with 
disk, envelope, and outflow cavities. The interstellar conditions are used to assess 
whether the star formation was spontaneous or induced by external pressure. 

Subject headings: Magellanic Clouds — HII regions — ISM: individual (N51) 

1. Introduction 

The Spitzer Space Telescope, with its high angular resolution and sensitivity in the 
mid- and far-infrared wavelengths, opens a new window to view the formation of individual 
massive stars not only in the Galaxy but also in nearby galaxies, such as the Large Magellanic 
Cloud (LMC). Using Spitzer observations, Jones et al. (2005) reported four young stellar 
objects (YSOs) in the H II complex N 159 in the LMC. We have obtained Spitzer observation 
of the superbubble LHa 120-N51D (designation from Henize 1956), or N51D for short, in the 
LMC, and are surprised to find YSOs projected within the superbubble where the density is 
expected to be low. 
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N51D has been partially imaged with the Hubble Space Telescope (HST) in the Ha 
and [S II] lines (Chen et al. 2000). These HST images allow us to identify dust globules 
and Herbig-Haro (HH) objects associated with YSOs. We have modeled the spectral energy 
distributions (SEDs) of these YSOs in dust globules to estimate their physical parameters. 
We have further used Ha and X-ray observations to determine the external gas pressure of 
the dust globules, in order to assess whether the star formation was spontaneous or induced. 
This Letter reports our analysis of Spitzer and HST observations of the three YSOs embedded 
in dust globules in N51D. 



2. Observations and Data Reduction 

Our Spitzer observations of N51D were made with the Infrared Array Camera (IRAC; 
Fazio et al. 2004) on 2004 December 16 and the Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004) on 2005 April 3. The IRAC observations were obtained using 
the mapping mode to cover a ~30'x20' area in the 3.6, 4.5, 5.8, and 8.0 /im bands. The 
four band maps overlap in a smaller region, ~23'x20' in size. Each position in the map 
was observed with five 30 s integrations taken in the high dynamic range mode and using 
medium dithers in a cyclic pattern between exposures to aid in the rejection of transients. 
The resulting integration time at each location in the image was ~150 s. Our analyses use 
the post-BCD (Basic Calibrated Data) products resulting from standard pipeline processing 
(version SI 1.0.2). 

The MIPS observations were obtained using the scan map mode in the 24, 70, and 160 
/im bands. The medium scan rate was used to map a region comprised of sixteen 0?5 scan legs 
with a cross-scan step of 148" to cover a region 30'x20' in all three MIPS bands. The MIPS 
DAT v3.00 (Gordon et al. 2005) was used to do the basic processing and final mosaicking of 
the individual images. In addition, extra processing steps on each image were carried out. 
At 24 /zm, these include readout offset correction and division by a scan-mirror-independent 
flat field. At 70 and 160 /im the extra step was a pixel-dependent background subtraction 
for each map (using a low-order polynomial fit to the source free regions). The final mosaics 
have exposure times of roughly 200, 80, and 18 s pixel -1 for 24, 70, and 160 /im, respectively. 
The MIPS photometry was determined using the point source fitting software StarFinder 
(Diolaiti et al. 2000) and a smoothed STinyTim model point spread function. 
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3. YSOs in the Superbubble N 51D 

3.1. Identification of YSOs 

YSOs can be distinguished from normal stars because of their excess infrared emission 
from circumstellar dust. To show the locations of YSOs and normal stars with respect to 
the ionized interstellar gas, we present in Figure la a color composite of IRAC 8 /im (red), 
optical continuum at 5000 A (blue), and Ha (green) images. In this color image, the ionized 
interstellar gas appears green and diffuse, most normal stars appear blue and unresolved, 
YSOs appear red and unresolved, and the interstellar dust emission appears red and diffuse. 

The ionized gas shell of the superbubble N51D is elongated along the N-S direction 
with a faint extension to the southwest. The 8 /im diffuse dust emission is the brightest 
in filaments along the periphery of the ionized gas shell; it also shows a faint component 
extending across the entire field. Two OB associations, LH51 and LH54 (Lucke & Hodge 
1970; Oey & Smedley 1998), exist within the superbubble, appearing as concentrations of 
blue stars near the western and eastern shell rims in Figure la. Bright YSOs are seen within 
the boundary of the N51D superbubble. Some YSOs along the north and west shell rims 
are associated with bright filamentary dust emission, while the other YSOs appear isolated. 
Some YSOs are projected within the OB associations LH51 and LH54, but the YSOs on 
the north rim are not near either association. 

Note however that not all red stars in Figure la are YSOs. These red stars may appear 
red or blue in the color composite of IRAC 3.6, 4.5, and 8.0 fj,m images displayed in Figure lb. 
To determine the nature of these stars, we examine their colors quantitatively and present 
an IRAC color-color diagram ([3.6] — [4.5] versus [5.8] — [8.0]) and a color-magnitude diagram 
([8.0] versus [3.6] — [8.0]) of N51D in Figure 2. Similar diagrams presented by Allen et al. 
(2004) and Whitney et al. (2004) for star forming regions show that normal stars have zero 
colors and that YSOs exhibit a range of red colors depending on their classes, or amounts of 
circumstellar material. For the color schemes we use, YSOs appear red in both Figures la 
and lb, but normal stars appear blue in Figure lb even if they appear red in Figure la. 

We use the objects 1-4 marked in Figures 2 and 3 as examples. Objects 1-3 are YSOs 
and they appear red in both Figures la and lb. Object 4 has been cataloged as L54S-81 
with V = 16.82±0.01, (B- V) = 1.67±0.05 (Oey 1996), and K = 11.64±0.02 (2MASS). Its 
colors and magnitudes are consistent with a Galactic M4 V star at a distance of 67 pc; thus 
star 4 appears blue in Figure lb and shows nearly zero colors in Figure 2a. 
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3.2. Dust Globules and Herbig-Haro Object 

The eastern rim of the N51D superbubble has been imaged by the HST WFPC2 in 
the Ha and [S II] AA6716, 6731 lines (Chen et al. 2000). Figure 3a shows a color composite 
of the HST Ha (green) and [S II] (blue) images and Spitzer 8 /im image (red). While the 
foreground M dwarf (object 4 in Fig. 2) appears isolated, the three YSOs (objects 1-3 in 
Fig. 2) are each associated with a small dust feature whose surface gas is ionized. 

Figure 3b shows a close-up Ha image of the dust globules associated with YSO-1 and 
YSO-2. The dust globules are each about 5" in size, corresponding to 1.25 pc, slightly larger 
than typical Bok globules in the Galaxy and in the LMC (Clemens et al. 1991; Garnett 
et al. 1999). The bright star to the east of YSO-2 is the WR star HD 36402 (Breysacher, 
Azzopardi, & Testor 1999). 

The dust feature associated with YSO-3, shown in Figure 3c, has a different morphology. 
Its surface toward the center of the N 51D superbubble is photoionized, but its back side 
is not. The [S II] close-up in Figure 3d shows additional knots of emission extending from 
the north rim outward along an almost straight line. These knots are faint in Ha. Their 
[S II] /Ha ratios increase outward from 0.5 to 0.65, significantly higher than those of the 
photoionized surface of the dust feature, 0.25±0.05. The alignment of these knots and their 
high [S II] /Ha ratios are characteristics that are frequently seen in HH objects (Reipurth & 
Bally 2001); thus we suggest that this is first HH object detected in the LMC. 

3.3. Physical Properties of the YSOs 

We have measured the flux densities of the three YSOs in N51D in all IRAC and 
MIPS bands, and listed them in Table 1. The SEDs were modeled with a 2-D radiation 
transfer code appropriate for YSOs (Whitney et al. 2003). These models include thermal 
dust emission but not PAH emission. The circumstellar geometry consists of any combination 
of a rotationally flattened envelope, flared disk, and bipolar cavities. YSO-1 is fit with a 
source luminosity of 10,500 L Q appropriate for a B2-B3 main sequence star, an envelope 
infall rate of 2 x 10" 4 M Q yr -1 , disk mass of 0.1 M , and bipolar outflow cavity opening 
angle of 15°. The disk and envelope centrifugal radius is 300 AU. These model parameters 
suggest a fairly young protostar. The disk radius could be varied by 50% and, with minor 
variations in other parameters, result in a similar SED. The total envelope mass is ~ 7OOM , 
but this is somewhat dependent on our choice of outer radius and the radial density power 
law. The disk mass is also uncertain since the long- wavelength contribution blends with the 
envelope. However, without a disk, the envelope would have a steeper density profile and 
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less near- and mid-IR flux density. Our best model SED for YSO-1 is shown in Figure 4. 
Ten inclinations are shown, and the best-fit inclination is ~ 50°. 

YSO-2 is fit by a similar model as YSO-1: an illuminating B3 star (2500 L ), an envelope 
mass of ~210 M Q and infall rate of 4 x 1O~ 4 M yr -1 , disk mass of 0.1 M and radius of 300 
AU, and bipolar outflow cavity opening angle of 30°. The best-fit inclination is ~ 60°. 

YSO-3 is interesting because it requires a massive envelope to fit the long-wavelength 
SED, but a nearly unextinguished line of sight to the illuminating central source, suggested 
by the optical data available for this source (see Fig. 4; Zaritsky et al. 2002). We fit this 
with a ~B4 star (1200 L Q ), envelope mass of ~230 M Q and infall rate of 5 x 1O~ 4 M yr -1 , 
disk mass of 0.2 M Q and radius 400 AU, and bipolar cavity opening angle of 40°. The 
best fit inclination is ~ 40°. This suggests either an evolved protostar with a substantial 
envelope and a large cavity carved out by outflows; or more than one source in the beam, 
an unextinguished source surrounded by a disk, and an embedded source with sufficient 
envelope mass to fit the long-wavelength SED. In either case, the model requires substantial 
warm dust, such as from a disk, to fit the mid-IR flux density. In the latter case, the HH 
object could come from either source. 

In all of these sources, PAH emission in the [3.6], [4.5], and [8.0] IRAC bands is likely to 
contribute additional flux densities over the thermal emission accounted for by the models. 
This is especially true in YSO-3 which is less embedded than the other two sources. We have 
thus excluded these bands in estimating the best-fit inclination. 

4. Spontaneous or Induced Star Formation? 

The three dust globules that harbor YSOs in N 51D are exposed to the radiation field 
of OB associations, similar to the dust pillars in the Eagle Nebula (Hester et al. 1996). The 
interstellar medium in N 51D has been well-studied, so we can use the external conditions of 
the dust globules to assess whether the star formation was induced by external pressure. 

The interior of a superbubble is expected to be filled with shock-heated gas at 10 6 -10 7 
K (Chu et al. 1995). The existence of hot gas in N51D has been confirmed by the diffuse 
X-ray emission detected by Einstein and ROSAT (Chu & Mac Low 1990; Dunne, Points, 
& Chu 2001). Recent XMM-Newton observations provide high-quality X-ray spectra, and 
the best-fit thermal plasma emission models indicate that N51D's interior hot gas has an 
electron temperature of ~ 2.6 x 10 6 K and an electron density of ~ 0.03 cm -3 (Cooper et 
al. 2004). The thermal pressure of this hot gas is P/k ~ 1.5 x 10 5 cm -3 K, assuming a He 
to H number ratio of 1:10. 
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The ionized gas enveloping the dust globule around YSO-1 has a peak Ha surface 
brightness of 1.4 x 10~ 14 ergs s~ x cm -2 arcsec -2 along its limb-brightened rim, corresponding 
to an emission measure of 6740 cm -6 pc for a temperature of 10 4 K. The ionized gas has a 
shell structure, and the longest path viewed through the shell is along the rim and equals 
2^ 2 R(AR/R) 1 / 2 , where R and AR are the shell radius and thickness, respectively. For 
a radius of 0.65 pc and a thickness of 0.15 pc (measured from the Ha image), the peak 
emission measure leads to an electron density of 90 ± 10 cm~ 3 . As the temperature of a 
photoionized gas is of order 10 4 K, the thermal pressure of the ionized skin of the dust 
globule is P/k ~ 2 x 10 6 cm~ 3 K, substantially higher than that of N 51D's interior hot gas. 

The thermal pressure in the dust globule is P/k ~ 10 4 cm -3 K, assuming typical values 
of density (10 3 H 2 cm -3 ) and temperature (10 K) for Bok globules. This pressure is much 
lower than those of the ionized surfaces or the surrounding hot gas. The radii of the dust 
globules are slightly larger than the Jeans radius of 0.5 pc, implying that spontaneous collapse 
may have started the star formation. However, their co-location with the ~ 3 Myr old OB 
association LH54 (Oey & Smedley 1998) suggests that star formation in the globules is 
delayed and may have been induced recently by the thermal pressure of the superbubble 
interior. Assuming additional magnetic support, the magnetic pressure in the dust globules 
would be <4 times the thermal pressure, suggesting a magnetic field of <35 /xG, compatible 
with those observed in Galactic molecular clouds (Crutcher 1999). 

5. Summary 

Using Spitzer and HST observations of the superbubble N51D in the LMC, we have 
identified three YSOs in dust globules. One of these is associated with an HH object that is 
detected outside the Galaxy for the first time. The SEDs of these YSOs are consistent with 
those of early-B stars with disk, envelope, and outflow cavities. The thermal pressure of the 
dust globule is much lower than that of the surrounding warm photoionized gas or the hot 
shock-heated gas in the superbubble interior. If star formation is induced by the thermal 
pressure of the hot gas, the magnetic field of the dust globule is <35 /iG. 

This work is supported by the grant JPL- 1264494. 
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Table 1. Spectral Energy Distribution of Three YSOs in N51D a 





Source 1 


Source 2 


Source 3 


Wavelength 


Flux 


Flux 


Flux 


(//m) 


(mJy) 


(mJy) 


(mJy) 


3.6 


4.16± 10% 


1.51± 10% 


2.14± 10% 


4.5 


9.35± 10% 


1.44± 10% 


2.35± 10% 


5.8 


24.7 ± 10% 


4.71± 10% 


5.35± 10% 


8.0 


47.6 ± 10% 


11.4 ± 10% 


13.3 ± 10% 


24 


356 ±10% 


43 ±10% 


20 ±10% 


70 


1508 ±20% 


429 ±20% 


180 ±20% 


160 


1402 ±20% 


<193 


<193 



a The errors are dominated by uncertainties in the flux cal- 
ibration, as given in the IRAC and MIPS Data Handbooks. 
The measurement errors are only 2-3% for sources 1 and 3, 
and 4-7% for source 2. 
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Fig. 1. — (a) A color composite of N51D with Ha in green, green continuum in blue, and 
IRAC 8 /im in red. The Ha and green continuum images are from the Magellanic Cloud 
Emission-Line Survey (Smith et al. 1999). (b) An IRAC color composite of N51D with 3.6 
/im in blue, 4.5 /zm in green, and 8.0 jum in red. 
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Fig. 2. — foj [3.6] — [4.5] versus [5.8] — [8.0] color-color diagram of N51D. (b) [8.0] versus 
[3.6] — [8.0] color-magnitude diagram of N51D. Objects 1-3 are marked by circles around the 
data points, and object 4 is marked by a star symbol. 



Fig. 3. — (a) A color composite of HST and Spitzer images of N 51D with Ha in green, [S II] 
in blue, and 8.0 /im in red. Note that the Spitzer 8 /xm image has a ~ 2" resolution so the 
sources appear as smudges, as opposed to sharp point sources in the HST images, (b) HST 



WFPC2 Ha image of YSO-1 and YSO-2. (c) HST WFPC2 Ha image of YSO-3. (d) HST 
WFPC2 [S II] image of YSO-3. 
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Fig. 4. — Models of SEDs for YSO-1 (left) and YSO-3 (right). The different colored lines 
correspond to different inclinations ranging from cost = 0.05 (green) to 0.95 (pink) in steps 
of 0.1. The observations are shown as triangles with error bars. 
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